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Abstract 



►^ ^ Einstein equations for several matter sources in Robertson- Walker and 

^^ • Bianchi I type metrics, are shown to reduce to a kind of second order nonlin- 

(N 

^ i ear ordinary differential equation y + af{y)y + (3f{y) J f{y) dy + jf{y) = 0. 

c/3 ' Also, it appears in the generalized statistical mechanics for the most inter- 
O 



esting value q = —1. The invariant form of this equation is imposed and 
the corresponding nonlocal transformation is obtained. The linearization of 



^ ■ that equation for any a, (3 and 7 is presented and for the important case 



/ = 6y" + k with /? = a"^ , ^+2\'^ i^s explicit general solution is found. More- 
over, the form invariance is applied to yield exact solutions of same other 
differential equations. 
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I. INTRODUCTION 

Exact solutions of the Einstein equations are difficult to obtain due to their nonlinear 
nature. There exist several interesting physical problems where the Einstein field equations 
for homogeneous, isotropic and spatially flat cosmological models with no cosmological con- 
stant [jl|- and for a time decaying cosmological constant [0, or Bianchi I-type metric |^ 
with a variety of matter sources, reduce to particular cases of the second order nonlinear 
ordinary differential equation 

y + af{y)y + (3f{y) j f{y) dy + ^f{y) = 0, (1) 

where y = y{x), f(y) is a real function and the dot means differentiation with respect to x. 
a, j3 and 7 are constant parameters. 

Recently, it was shown that some galactic models of astrophysical relevance, when investi- 
gated with the "generahzed" Statistical Mechanics [0] , can be exactly described by solutions 



to the Boltzmann equations that maximize the Generalized Tsallis Entropy for g = — 1 [10 



and it was found that the corresponding probability distribution function satisfles (|I]) |11[ |. 
It is believed that quantum effects played a fundamental role in the early Universe. For 
instance, vacuum polarization and particle production arise from a quantum description of 
matter. It is known that both of them can be modeled in terms of a classical bulk viscosity 



|T^. Using the relativistic second-order theory of non-equilibrium thermodynamics called 



Extended Irreversible Thermodynamics developed in |T^ [|^, it was considered a homoge- 
neous isotropic spatially-flat universe, flUed with a causal viscous fluid whose equilibrium 
pressure obeys a 7-law equation of state, while the transport equation of the viscous pressure 



IS 



a + Ta = -SCH-^era Uh + ^-^-^V (2) 



with e = [|l^. Following ||T6[ for m = 1/2, it was shown in [|l| that the expansion rate 



satisfies a modified Painleve-Ince equation that has the form of (||) with f{y) = y and 7 = 0. 
Cosmological models with a viscous fiuid source have been studied using the full causal 
irreversible thermodynamics with the full version of the transport equation for the bulk 
viscous pressure [|T^ P 0. Relating the equihbrium temperature T with the energy density 
in the simplest way to guarantee a positive heat capacity, it was shown that the expansion 
rate satisfies (|T]) for m = 1/2, with f{y) = y~'^^^ and 7 = 00. Also, the early time evolution 
of a dissipative universe, leads to an equation for the expansion rate that has the form (|l]) 
0], [|18|, in the relaxation dominated regime. 



Another interesting example appears when an anisotropic universe, described by a 
Bianchi type I metric, is driven by a minimally coupled scalar field with an exponential 
potential. The Klein-Gordon equation for the scalar field and the Einstein equations for 



the metric are expressed in term of the semiconformal factor G and their derivatives p9 
Then, the solutions of this equation set can be obtained if one is able to solve the following 
Einstein equation for G, 

g| + (c-1)G+| = C2, (3) 

which, making the substitution G = y^^^ (^ becomes (|I]) [p. A similar result is obtained 
in the particular case when the Bianchi type I metric reduces to a fiat Robertson- Walker 
space-time 0. 

From the Generalized Tsallis Entropy, defined as ^ 

^. = %-i)-^Efe-pD, (4) 



it can be constructed the generalized Statistical Mechanics where fc is a positive constant, q 
is a real number that characterizes the statistic and the sum is made over all the microscopic 
configurations whose probabilities are pi. It leads to the conventional Boltzmann-Shannon 
statistic in the limit g — > 1 and it is found to be a good framework to study astrophysical 
problems, as the Generalized Freeman Disk [^ and Kalnajs oscillations of a slab of stars 
| 2T| . Taking the generalized Fisher information for Tsallis Statistics [^ 



— f 



2 



where /d(x) is the probability distribution function, and solving the variational problem in 
order to find the distribution function that maximizes the Fisher information, a differential 
equation of type (|I]) is obtained for y = fd/fd, where f{y) = y, a = (2g — 1), /5 = |g(g — 1) 



and 7 = ||rT| . For relevant physical applications the most interesting value of the statistic 



parameter is g = — 1 |]T0|, in this case the above equations can be solved explicit and the 
general solution will be given in section 3. 

Thus, it turns out to be of great interest to analyze ([^) from the physical and mathemat- 
ical point of view. The paper is organized as follows, in section II we introduce an invariant 
form and use it to reduce (|l|) to a linear, inhomogeneous ordinary second order differential 
equation with constant coefficients, by means of a nonlocal transformation. Then, its para- 
metric general solution is given. In section III we extend the nonlocal transformation and 
find the explicit general solution of a modified Painleve-Ince equation for /5 = 1/9 I^Sj. In 
section IV we use the nonlocal invariance to obtain a new class of differential equations for 
which the general solution is found. In section V the conclusions are stated. 



II. FORM INVARIANCE 



The differential equation (|ID, which appears in several interesting physical problems, 
has been solved and studied in particular cases using nonlocal transformations as it was 
previously stated. To investigate (|I]) we write it in invariant form 



+ ay + (3 j f{y)dy + ^ = j— + ay' + (3 j f{y)dy + -i, (6) 

under the nonlocal transformation group defined by the transformation 

Pf{y)dy = p7{y)dy, (7) 



-fiy)dx = -f{y)dx, (8) 

a a 



£ = £. 



(9) 



/9c + 7 = /3c + 7, (10) 

where f{y) is a real function of y = yix)-, the prime indicates differentiation with respect 
to X. a, /3, 7 are constant parameters and cic) is an integration constant provided by the 
integral on the left (right) hand side of (P). By invariant form we mean that the left hand 
side of (ID transforms into the right hand side under the nonlocal transformation defined by 
(0^10) f*^^ ^^y functions /, /. The parameters a, /3, 7, a and [3 satisfy (|9|-p^. 

The form invariance group can be used to linearize (|l]). In fact, taking the function 
fijj) = 1, a = a, P = P and 7 = 7 (this means c = c) in the invariant form (|) and the 
transformation ([7|-p!0D, they become 



fiy) 



+ ay + [3 j fiy) dy + j = y" + ay' + py + pc + -f, (11) 



y = jfiy)dy, x = jfiy)dx. (12) 

Without loss of generality we choose c = c = 0. So, if the invariant ( pT]) vanishes, then, (|l]) 
transforms into 

y" + ay' + Py + ^ = 0, (13) 

under the transformation of variables (0). This is a linear, second order ordinary differential 
equation with constant coefficients . Its general solution is 



a)/3^ 



y = ci exp (Aix) + Ca exp (Aax) - -, (14) 



where Ai and A2 are the roots of the characteristic polynomial of (|T3D. We indicate the 
integration constants with c, Ci,...c„ and c, ci,...Cn- 



2 



y={ci + C2x)expi^--j --. (15) 

The real solutions can be classified as follows (we also assume that a, (3 and 7 are real). For 

2 2 

a > and /3 < ^ we have two real, negative roots for a strong damped solution. For (3 = ^ 

2 

we have a double-negative root for a critically damped solution. For a > and /? > ^ we 
have two complex roots with negative real parts for a weakly damped solution. For the case 
a < growing solutions occur. 

The transformation of variables (p!2|) , relates the general solution of @) with y{x) through 



(|1|). We find that 



X 



y = y{y{x)), (16) 



■dx. (17) 



fiyivim 

are the parametric equations for x and y in terms of x. In the particular case f{y) = y we 
have shown that a class of nonlinear modified Painleve-Ince equation can be transformed 
into a linear second order ordinary differential equation by a nonlocal transformation. 

The theory introduced by Lie considers the invariance of the differential equations under 
point transformations. He showed that the one- dimensional free particle equation has the 
eight-dimensional SL(3,R) group of point transformations. This is the maximum number of 
symmetry generators for a second-order differential equation of the form |^ 



y + h{y,y,x) = 0. (18) 

In our case (|lD has the form of (0). Then, it has eight or less point symmetries. However, it 



becomes (13) under the transformation of variables (12) and can be cast into the free particle 
equation by a local point transformation. So, ([T3|) has always eight symmetry generators. 
We conclude this section observing that the nonlocal transformation ([?|-|T0|) changes the 
number of symmetry generators for the class of differential equations (|l]) and the physics 
contained in the original problem. 

A. The nonconstant parameters case 



Here we allow the parameters in (|T|) and in the transformation ([7|-|T0|) to be functions of 
the independent variable, that is, a = a{x), (3 = l3{x) and 7 = 7(x). In order to preserve 
the form (|I|) we choose a{x) = a{x) and P{x) = P{x). In this case, the invariant form (^ 
reads 

8 



-^ + a{x)y + /?(x) j f{y) dy + 7(x) = 4=7 + «(3^)F' + /5(3^) / fiV) dy + l{x), (19) 

where x is the transformed of the point x. Therefore, taking 7 = 7 and f{y) = 1 we can 
hnearize the equation 



y + a{x)f{y)y + I3{x)f{y) / f{y) dy + i{x)f{y) = 



(20) 



which transforms into 



y" + a{x)y' + [3{x)y + 7(x) = 0. 



(21) 



An important physical problem of general relativity, concerning the motion of expanding 
shear-free perfect fluids [^, is governed by the ordinary differential equation 



y = F{x)y'^ 



(22) 



where F{x) is an arbitrary function from which the equation of state can be computed. A 
complete symmetry analysis of this differential equation was given in |2B]. Here we see that 
it is contained in the set of equations ( PP] ) when a{x) = 0, /3{x) = ^^^ , 7(0;) = and 
f{y) = y^^'^ ■ Then, choosing /(y) = (y)^^/^ in (|7|- |lO|) , the transformation of variables is 



y = — , X = — dx, 
^9 7 3 



(23) 



and (|2^) becomes 



y" = 3F(x), 



(24) 



thus 



Fix) dx 



dx + cix + C2, 



(25) 



is the general solution of the simple linear equation (p 



III. EXTENDED NONLOCAL TRANSFORMATION 



The integral in (|T^ can be performed analytically and the general solution y = y{x) of 
(H) obtained explicitly for a special set of functions f{y). For this purpose we generalize the 
nonlocal transformation group defined by ([^-|T0|) extending it to 



fuiy) dy + fi2{y) dx = fuiy)dy+ fi2iy)dx, 



(26) 



f2i{y) dy + f22{y) dx = f2i{y)dy+ f22{y)dx. 



(27) 



For simplicity we begin our investigations restricting ourselves to the case x = x, that is, 
/21 = /21 = 0, /22 = /22 = 1 ^'^d requiring the invariant form (H) to be invariant under the 
remaining nonlocal transformation group, defined by (^6]-^) with the above restrictions. 
Under these assumptions we can write the nonlocal transformation as 



y = p + qy, 



(28) 



where the functions p and q are expressed in terms of the functions /n, /12, fn and f^2- So, 
they have a specific dependence on the variables y and y 



p{y,y) 



fi2{y) fui 
7u{y) 7ii( 



Qiy,y) 



hi{y) 

fniV) 



Inserting (^Hj) in (||) we get 



(29) 



(30) 



+ ay + pjfdy + -f = =y + 



dq dq 
dy dy 



y' 
/ 
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dp dp dq 
oy oy oy 



y p dp _. ., -n I -p ,_ _ 

f fdy 



(31) 



and comparing the coefficients of y , we have 



dq dq 
oy dy 



(32) 



whose solution is 



(i{y,y) 



(33) 



Using (|33D and comparing the coefficients of ij we easily find that f = y and f = y. But, 
the comparisons of the coefficients of y and the remaining terms give the equations 



a 



dp y dp p 
dy ydy y 



1 _y 

- + a-, 

y y 



(34) 



p dp 
/3 jydy + ^ = --^ + ap + [3 / ydy + 'y. 



(35) 



The function p that satisfies ( p4|) is given by 



(y. Oi 

p{y, y) = -^yy - ^v^ + h{y^y)^ 



(36) 



where the function h{y, y) satisfies the partial differential equation 



dh dh ^ ^ 
dy dy 



(37) 



It can be seen that the solutions of (pTf ) are given hy h = ho/y, where ho is an arbitrary 
function of the quotient y/y. So, the form of the solution for p is 



. _N " _ "_2 , ho{y/y) 
p{y,y) = ^yy- ^y + • 



(38) 



Comparing (|30D with (|33| ) we have fniy) = 1/y, and comparing (|29|) with (^), we obtain 
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ho{y/y) = C1- + C2-. 
y y 



Inserting ( p9D in (|35D we find that ci = C2 = 0, 7 + /3c = /5 c + 7, and 



(39) 



P 



2a' 



'7=^2 



P 



2a 



(40) 



Therefore, the final invariant form and the resulting nonlocal transformation are 



y 
y 



a 
"9 



9 



^ . " 2 /I y —- Ct _2 -n- — 

- + ay + -—y + IJC + 'j = - + ay + -—y +pc + 7, 



y a y a_ 
1/3 y 3 



(41) 



(42) 



In the particular case in which the invariant form (^) vanishes, the l.h.s. gives rise to a 
nonlinear differential equation 



a 



y + ayy + yi/^ + 72/ = o. 



(43) 



(where, without loss of generality we have taken c = c = 0, so that, 7 = 7), that can be 
solved using the invariance properties formulated above. To do this, we make a = on the 



r.h.s. of (^ll) . Then, inserting its solution in (^21), it can be integrated giving the general 
solution 



y 



3 2cix + C2 
a cix^ + C2X + C3 ' 



7 = 0. 



3 /^ ci exp (^/^x) + C2 exp (-y/ix) 
y=— / N r^ v^ — , 7 7^0. 



(44) 



(45) 



" ci exp (v^x) - C2 exp [-^/yx^ + C3 ' 
It can be seen that (|i3| ) has eight Lie point symmetries and it is equivalent to a second order 



linear differential equation under a point transformation |^^. On the other hand, for any 



other value of the coefficient j3 7^ 



2al 
9 ' 



has two point Lie symmetries and we cannot find 
12 



a point transformation that cast it in a linear equation ||2^. However, using the invariant 
form ( |Tl| ) and the transformation of variables (|T2|) for f = y, we have proved that ( ^3[ ) 
can always be linearized whatever the value of the coefficient of y^ is. Therefore, using the 
invariance properties of the form (^ we have obtained the same results that come by the 
Lie theory of symmetries. In addition, we have linearized (HI) when it has less than eight 



Lie point symmetries. 



IV. SOLUTION OF NEW CLASSES OF DIFFERENTIAL EQUATIONS 



Now, we are going to investigate the case when the invariant expression (^) vanishes, 
and we shall construct several important classes of solvable second order nonlinear ordinary 
differential equations. To do this, we must seek the nonlocal transformation defined by 
(^,^) with the condition that the invariant (^Tj) vanishes. This leads to the equations that 
determine it 



r y 
af = - 

y 



dp y dp p 
dy ydy y 



+ «/, 



(46) 



Pf I fdy + if = l 



P^ + apf + P'f I fdy + -ff 



(47) 



and we shall show a set of functions /, / for which the nonlocal transformation exists. The 
solution of (|46|) can be obtained writing 



p{.y,y) = aypo{y){y) + pi{y) +P2{y,y), 



(48) 



where each function satisfies 



f = 2po + ypQ, 
13 



(49) 



Pi + — + af = 0, 

y 



dp2 , _dp2 , „ 



(50) 



(51) 



where the ' indicates derivative with respect to the argument of the function. Solving the 
system (i^-^) and inserting their solutions in (^Sj), we find the solution of (^), that is: 



Piy,y) = tt^ yfdy- - / fydy + — . 

y J y J y 



(52) 



Comparing (^2|) with (P^D, the function ho{y/y) is given by (^), but these terms can be 
absorbed in a redefinition of the integration constants provided by the two integrals of 



Then, without loss of generality we take them equal to zero. 

From ( ^|52D we obtain the difficult integrodifferential equation that satisfy the functions 
/ and /. It reads 



a 






I fydyV + (3^- I fdy + ^L 



fydy 



f 



, / 



+ /5= fdy + j^. 

y J y 



(53) 



In what follows we shall show a set of functions /, / that are solutions of this integrodif- 
ferential equation and construct three sets of nonlinear differential equations that can be 
linearized and explicitly solved. 

A. Case a 



An interesting solvable equation set can be obtained when we choose the functions /, / as: 



/ = %" + A;, f = hy^ + k. 



(54) 
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Taking into account that the left hand side of (|53| ) depends of y and its right hand side 
depends of y, it must be a constant. So, inserting the functions given by ( |5^ ) in (^) and 



after some algebra, it provides the constrains satisfied by the parameters 



(3 = a 



2 ^+ 1 

(n + 2)' 



(3 = a 



2 n + l 



(55) 



/3A;2-a2^=^fc2_^2^^ 

4 4 



In addition, the function piy, y) is given by 



p{y, y) = ay 



b n k 

y +- 

n + 2^ 2 



ay 



n 



-r + 2 



Finally inserting ( ^^ - |56|) in the invariant form (|^), we have 



i/ + « [6y" + fc] ?/ + /3 



uiV 



2n+l 



6A; 



n + 2 



n + 1 n + 1 



l/'^^' + ry 



(56) 



(57) 



(58) 



^ + a by'^ + k ij + (3 



py'''^' ,7:t:^+2...+i 



n+ 1 



+ 6A; 



n + 1 



y 



+ /c y 



(59) 



Besides, from ( p8| , |33| , |57D we obtain the nonlocal transformation (^) in invariant form 



V aby^ ak v aby^ ak 

- H \ = - H \ 

y n+2 2 y n+2 2 ' 



(60) 



that links (|58|) and (^). To integrate these equations we use their invariant property along 
with (|55| - |56D and analyze two different cases. In the first case, we choose b = 0, a = a, 
k = k and n = n. Then, (3 = (3 hj (^) and (|59|) reduces to a linear second order differential 
equation ioi y = y with constant coefficients 



" 7 - 2 7 2 '^ + 1 

y + aky + a k - — — — 
(n + 2Y 



y = 0. 



(61) 



Integrating (^) for the above value of the parameter, we obtain the general solution of (^ 
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n + 2 y" 
abn J y"- dx^ 



(62) 



where y is any solution of (|6TD . In the second case, when we choose b = 0, a = a, k = k 
and n = n, the (|58|) reduces to (|6lD ior y = y and the general solution of (|59| ) is 



_^_n + 2 y^ 



(63) 



a 6 n J y"- dx 

where y is any other solution of (|6T|). Inserting the general solution of the (|58D and (p9D, 
given by ( |62D and (|63D, in the nonlocal transformation (pO|), it can be integrated and the 



final relation between the variables y and y, that transforms ( |58| - p9D one on each other, is 



y"' dx 



ak ^ 
exp -X 



y"- dx 



(ak ^ 
exp — X 



(64) 



For the particular case n = n = — 1, we obtain 7 = ah and 7 = ah. All the remaining 



equations (|60D-(|64D can be applied for ra = — 1 and n = —1 because they do not depend 
explicitly of the parameters /?, /3, 7 and 7. 



In the next subsections we investigate other generalizations of (|58| - |59D , that can be lin- 
earized and solved. 



B. Case b 



Writing the equations set ( ^81) and (]59|) as 



F{y,y,y) = 0, F{y,y,y) = 0, 



(65) 



a generalization of both equations can be done expressing them in the following way. 



1 1— ^ ^ _ 

-F{y,y,y) = -F{y,y,y), 



(66) 
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which is invariant under the nonlocal transformation given by ([60|). It is easy to prove that 
the new functions 



F{y, y, y) = F{y, y, y) + 5y, F{y, y, y) = F{y, y, y) + 



(67) 



where 5 is a constant parameter, also satisfy the invariant condition 



1 ~ 1— ^ ^ _ 

-Fiy,y,y) = -F{y,y,y). 

y y 



(68) 



This gauge symmetry generates a new nonlinear equation that can be linearized and solved. 
In fact, when the invariant in (|B^) vanishes, it gives rise to a set of equations that transform 
one on each other under the same nonlocal transformation, these are: 



y + a [by'' + k]y + p 



aV 



2n+l 



n 



- + bk- 

n + 1 n + 1 



y"+^ + k'y 



+ 5y = 0, 



y + a 



by" + k t + /3 



b r + b k _ . ^ 1/"+^ + k y 



n + 1 



n + V 



(69) 



(70) 



In particular, to solve (^) we choose b = 0, a = a, k = k and n = n {f3 = (3 hj (|56|) ) in 



(ffOp. Then, it reduces to 



y + aky + 



2,2 ^ + 1 

a k 



in + 2y 



+ 5 



y = 0. 



(71) 



Inserting the solutions of (|7T|) in (|60| ) and integrating it for the selected parameters, we 
reduce (^) to quadratures 

'n + 2 y" 



abn Jy"'dx 



(72) 



For the particular case b = b = 1, k = k = 0, n = n = l and 6 = 'j, (p9| - [70D reduce to 
(^), the variable transformation ( pOD reduces to (^) and ( p^ ) gives the relation between 
the variables ?/ and y that leaves invariant (^l]) . 
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C. Case c 



There is an important result that can be deduced from ( |60D when a = a and k = k^ in 



this case the nonlocal transformation 



y_ aby 



is fc-independent, 

y aby" 



" 777 



y n + 2 y n + 2' 



(73) 



and by (p5|-|56|) 



n = n, 



-n 



n 



n + 1 



(74) 



So, if we take k{x) and 5{x) as functions of the independent variable x instead of constant 
parameters, then, there is no change in the deduction of the variable transformation ([73|), 
that comes from (^6|-^). This means that the set of equations (|69| - |70D give rise to new 
solvable equations that transforms between them by the nonlocal transformation (|73| ) 



y + a [by'' + k{x)] y + l3 



y + a by'^ + k{x) y + [3 



b^y _ + bk{x)^^-^y^+' + k\x)y 



n + 1 



n+V 



+ 5(s)y = 0. 



(75) 



b'y^ + bk{x)'^y''^' + k\x)y 
n+1 n+1 



+ 5{x)y = 0. (76) 



For instance, to obtain the solutions of ( |7SD we take 6 = and n = n in (|T6D and it becomes 
a general homogeneous linear second order differential equation 



y + ak{x)y + 



a k [x] 



+ 6{x) 



y = o, 



(n + 2)2 
then, inserting the solutions of this equation in (|7^) , we reduce (^) to quadratures. 



(77) 



V. CONCLUSIONS 



We have introduced a new invariance concept that leads to classes of second order non- 
linear ordinary differential equations which are equivalent under nonlocal transformations. 



These classes contain a second order linear ordinary differential equation with constant co- 
efficients. The parametric expression of the solutions for an arbitrary function f{y) and any 
values of the parameters a, (3 and 7, has been found. Also, the case in which these pa- 
rameters are functions of the independent variable has been investigated. Several important 
physical problems are mathematically described by these equation classes. Many of these, 
arise in General Relativity when the Einstein field equations are investigated for homoge- 
neous, isotropic and spatially fiat cosmological models with no cosmological constant, or 
Bianchi I-type metric with a variety of matter sources. Also, the probability distribution 
function, which maximize the Fisher's information measure in the generalized Statistical 



Mechanics, was found to satisfy (^) for the most interesting value g = — 1 ||TT |. 

Taking x = x m the nonlocal transformation, and imposing the form invariance of the 
general expression (^), we have obtained a modified Painleve-Ince equation (|43|). The non- 
local transformation of variables and the general solution of these equations has been found. 
In this case the equation has the eight dimensional group of Lie point group symmetries 
SL(3,R) and this is the maximum number of point symmetries that a second order differen- 
tial equation can have. Other sets of new nonlinear second order differential equations are 
generated, that can be linearized and solved explicitly (^,^,^). It is also to be remarked 
that, the use and application of the form invariance have lead to exact solution of differential 
equations whose solution were unknown, in particular for modified Painleve-Ince equations 
and polinomical differential equations, which usually appear in problem related with quan- 
tum effects in the very early Universe, originated by the vacuum polarization terms and 
particle production arising from a quantum description of matter, or when both of them are 
modeled in terms of a classical bulk viscosity 
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In general, the problem of finding solutions of nonlinear ordinary differential equations 
remains open. One direction along which one can proceed is to reduce them to a linear 
ordinary differential equation. For instance, when ([1|) possesses eight-parameter Lie group 
it is linearizable by a point transformation. On the other hand, the nonlocal transformation 
(0-[Ty) linearizes ([^) even when it has less symmetries. Thus, it could mean that has more 
nonlocal symmetries. We conclude that it is very interesting to study this kind of nonlocal 
transformations of variables and their associated nonlocal symmetries, which have received 
up to now little attention. We shall continue exploring this subject in future papers. 
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